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Abstract 

Using the non-equihbrium quantum field theory, photon production from 
the coherently oscillating axion field in a flat Robertson- Walker cosmology 
^ , is re-examined. First neglecting the Debye screening of the baryon plasma 

\ to photons , we find that the axions will dissipate into photons via spinodal 

I instability in addition to parametric resonance. As a result of the pseudo- 

■ scalar nature of the axion-photon coupling , we observe a circular polarization 

O . asymmetry in the produced photons. However, these effects are suppressed 

, to an insignificant level in the expanding universe. We then briefiy discuss 

a systematic way of including the plasma effect which can further suppress 
Q^' the photon production. We note that the formalism of the problem can be 

p • applied to any pseudo-scalar field coupled to photon in a thermal background 

1^ ■ in a general curved spacetime. 
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I. INTRODUCTION 



It is compelling that most of the matter in the universe is in a form of non-baryonic 
cold dark matter. If it exists, it would play an important role in the structure formation 
of the universe [|l|. Axions, the pseudo-Goldstone bosons, are among the most promising 
candidates for the non-baryonic cold dark matter. They arise from the spontaneous breaking 
of a global U{1) symmetry of Peccei and Quinn (PQ), which is introduced to solve the strong 
CP problem of QCD 0-^. In the standard big-bang cosmology, after the spontaneous 
breakdown of the PQ symmetry, the expectation value of the axion field (i.e. the axionic 
condensate) takes some random value on the interval [0, 2tt] and, is approximately constant 
over length scales which are smaller than the horizon size . If inflation occurs either after 
or during the PQ symmetry breaking, then the expectation value can be nearly constant 
throughout the entire universe 0. At high temperatures above the Aqcd scale, the axion is 
massless; however, at low temperatures, the axion develops a mass due to QCD instanton 
effects d^. Once the axion mass becomes greater than the universe expansion rate, the 
expectation value of the axion field begins to oscillate coherently around the minimum of its 
effective potential that is near the origin. The oscillating axion field then dissipates mainly 
due to the universe expansion as well as particle production 

In the original papers 0, simple estimates of the thermal dissipation of the homogeneous 
axionic condensate were given. They considered instabilities arising from the parametric 
amplification of quantum fluctuations that could pump the energy of the homogeneous 
axionic condensate into its quantum fluctuations via self couplings, as well as into quantum 
fluctuating photon modes via a coupling of the axion to electromagnetism due to the color 
anomaly of the PQ symmetry. This dissipational dynamics via quantum particle production 
exhibits the feature of unstable bands, and an exponential growth of the quantum fluctuating 
modes that are characteristics of parametric resonance. The growth of the modes in the 
unstable bands translates into profuse particle production. A given unstable mode will grow 
as long as it lies within the unstable band. However, eventually it will be red-shifted out 
of the band as the universe expands, and then the instabilities of parametric resonance 
are shut off. In Ref. 0, it has been shown that for the PQ symmetry breaking scale 
fa > lO^^GeV, because the axion is very weakly coupled, the time it takes to be red-shifted 
out of the unstable band is too short to build up an appreciable growth of the quantum 
fluctuating modes. Thus, all of these effects are insignificant. The condensate is effectively 
nondissipative and pressureless. It would survive in the expanding universe, and it behaves 
like cold dust at the present time. Interestingly, if fa ~ lO^^GeV, it could constitute a major 
component of the dark matter of the universe. 

Recently, the authors of Ref. |^ were motivated by the recent understanding of the im- 
portant role of the spinodal instability and parametric resonance that provide the nonlinear 
and nonperturbative mechanisms in the quantum particle production driven by the large 
amplitude oscillations of the coherent field |9Hl3|. They re-examined the issue of the dis- 
sipation of the axion field resulting from the production of its quantum fluctuations. They 
confirmed that the presence of the parametric resonance would lead to an explosive growth 
of quantum fluctuations if the universe was Minkowskian. Taking account of the expansion 
of the universe, quantum fluctuations of the axion do not become significant. This result 
confirms the conventional wisdom. 
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In this paper, we will re-examine the damping dynamics of the axion arising from pho- 
ton production in an expanding universe in the context of the non-equilibrium quantum 
field theory. The goal of this study is to present a detailed and systematical study of the 
above-mentioned problem using a fully non-equilibrium formalism ||9|-[13|. We will derive the 
coupled nonperturbative equation for the axion field and the mode equations for the photon 
field in a flat Robertson- Walker spacetime within the nonperturbative Hartree approxima- 
tion that is implemented to consistently take the back reaction effects into account. We then 
try to study both numerically and analytically how the nonperturbative effects of spinodal 
instability and parametric amplification of quantum fluctuations trigger photon production 
from the oscillations of the axion field. At this stage, it is worthwhile to mention that our 
approach can be generalized to any pseudo-scalar field coupled to the photon field in a more 
general curved spacetime. Since the pseudo-scalar nature of the coupling between the axion 
and the photon, the axion field affects the left- and right-handed circularly polarized photons 
differently. This leads to producing the two polarized photons in different amounts. This 
polarization asymmetry, if it survives, may have interesting effects on the polarization of the 
cosmic microwave background. 

To consider the fermionic plasma effect on photon production, one must systematically 
obtain the non-equilibrium in-medium photon propagators and the off-equilibrium effective 
vertices between the axion and the photon by integrating out the fermionic field to deal 



with this problem ||T2[. In a plasma, the transverse photons are dynamically screened [|14 
However, in the literatures 0, the arguments stated to include the fermionic plasma effect 
in support of their conclusions amount to adding by hand the electron plasma frequency 
into the propagating photon mode equations. This is problematic when we consider prop- 
agating photon modes in the presence of a thermal background. In fact, the consequence 
of the Abelian ward identities reveals that the transverse photons have the vanishing static 
magnetic mass in all orders of the perturbation theory [jl4|. This means that the in- medium 
transverse photon propagators must be nonlocal in nature, and cannot be approximated 
by the local propagator as suggested in Ref. even in the low energy limit [jl2|. Besides, 
in a fermionic plasma, the effective coupling between the axion and the photon resulting 
from integrating out the fermionic thermal loop can be modified both at finite temperature 
as well as out of equilibrium [|12|. Therefore, to fully consider the plasma effect, the non- 
equilibrium in-medium photon propagators as well as the off-equilibrium effective vertices 
play the essential roles. However, incorporating these non-equilibrium effects is a challenging 
task that lies beyond the scope of this paper, but certainly deserves to be taken up in the 
near future. In the following, we will totally ignore the fermionic plasma effect in order to 
focus on understanding whether the particle production due to spinodal instability as well 
as parametric amplification is effective or not in the cosmological context. 

In Sec. m, we introduce the axion physics and its coupling to the photon. An effective 
action of the axion-photon system in the expanding universe is derived. Sec. |ITT| is devoted 
to the formalism of the problem in terms of non-equilibrium quantum field theory. We 
obtain the equation of motion for the classical axion field and the photon mode equations. 
In Sec. |V[ we present the numerical results. Sec. is our conclusions. 
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II. EFFECTIVE ACTION OF AXION-PHOTON COUPLING IN EXPANDING 

UNIVERSE 



The physics of axion and its imphcation to astrophysics and cosmology can be found in 
the review articles 0,0. The axion does not couple directly to the photon. However, at 



tree level, the axion field has a coupling with the fermionic field ip, 

L<f,^ = g(j)ipi^. (2.1) 

Therefore, the axion can couple to the photon via a fermionic loop. As a consequence of the 
color anomaly of the PQ current, similar to the pion-photon system, the effective Lagrangian 
density for the axion-photon coupling is 

L^A = c^e'^^'^T^pF^,, (2.2) 

Ja 

where = df,A^ - d^A^,. In Eq. (|^), the scale fa = fpq/N, where fpg is the PQ 
symmetry breaking scale, and N is the color anomaly of the PQ symmetry. The coupling 
constant c = a^Epq/N — 1.95)/(167r), where Epq is the electromagnetic anomaly of the PQ 
symmetry, and a is the fine structure constant. Henceforth, we assume an axion incorporated 



into the simplest GUTs with Epg/N = 8/3, such that c ~ 1.04 x 10""^ p| 



Now we write down the effective action for the axion-photon system in the expanding 

S = J d^x^ (^L^ + La + -^L^a^ ' (2-3) 

where l/^/g is added in front of L^a, given by Eq. ( |2.2|) , because e°^/^^ is a tensor density 
of weight —1 ll^. In the following, for simplicity, we will assume a flat Robertson- Walker 
metric 

ds'^ = -g^^dx^'dx" = dt^ - a?{t)dx^, (2.4) 
where the signature is ( — h ++), and a{t) is the cosmic scale factor. In Eq. (|2.3|) , 

= -^g^-^d.^d,^ - V{(P,T), (2.5) 

La = -^g''>'g'"'F^f,F^,, (2.6) 

where the axion potential has a temperature-dependent mass term due to QCD instanton 
effects, being of the form 

V{<j>,T) = ml{T)f^(^l-cos^^, (2.7) 

m,(T)^0.1m„o(^)''', (2.8) 

where iriao is the zero-temperature axion mass, satisfying iriaofa — 6.2 x 10~^GeV^. Also, 
we use Aqcd - 200MeV. 
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It is well known that the minimal coupling of photons to the metric background is 



conformally invariant |T^. As such, in the conformally flat metric (|2.4| ), it is convenient to 
work with the conformal time, df] = a~^{t)dt. Hence, defining = x/a, the action ( p.3| ) 
becomes 



S= drid^x C 



Jdv 



d^x 



2 [dr]^ 



2 [dx 



+ - — —X -a V{-,T) 



2a drf 



4 a/a 



(2.9) 



where rj^^^ is the Minkowski metric. In terms of the conformal time, the effective action now 
has analogy with the effective action in Minkowski spacetime with the time dependent mass 
term and interactions. 



III. EQUATIONS OF MOTION 



The non-equilibrium effective Lagrangian in the closed time path formalism ||9HTT| is 
given by 

Cne, = C A+] -C[x-, A^] , (3.1) 

where + (-) denotes the forward (backward) time branches. We then decompose iiito 
the axionic mean field and the associated quantum fluctuating fields: 

X^{x,ri)=^{ri)+i)^{x,ri), (3.2) 

with the tadpole conditions, 

(7/>±(f,77)) = 0. (3.3) 

We will implement the tadpole conditions to all orders in the corresponding expansion to 
obtain the non-equilibrium equations of motion. 

To take account of the back reaction effects on the dynamics of the axion field from 
quantum fluctuating photon modes, we adopt the following Hartree factorization which is 
implemented for both ± components P-[TT[|: 



i^FF ipiFF). (3.4) 

As seen later, the expectation value can be determined self-consistently. It must be noted 
that there is no a priori justification for such a factorization. However, this approxima- 
tion provides a nonperturbative framework that allows us to treat photon fluctuations self- 
consistently [|ll| . On the contrary, we will ignore the quantum fluctuations of the axion. 



which can be produced via self-couplings, as the study of Ref. |^ has shown that these effects 
are insignificant. 

With Eq. ( |3.2D , we first expand the non-equilibrium Lagrangian density ( |3.1| ) in powers 
of ip and keep the term up to linear ip to ignore its quantum fluctuation effects. Together 
with Eq. (p.4|), the Hartree-factorized Lagrangian then becomes 
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a/a 



where 



a{ri) 



(p{ri) + a^{ri) m^(T)/a sin 



(3.5) 



(3.6) 



The dot means the time derivative with respect to the conformal time. 

With the tadpole conditions ( p.3|) , we obtain the following equation of motion for the 
axionic mean field given by 



a{r]) a^{ri) 



(3.7) 



where we define the dimensionless field, 9{rj) = ip[rj) / {a{vi)fa). 

Within the Hartree approximation, the photon production processes do not involve pho- 
tons in the intermediate states |11,13|. To avoid the gauge ambiguities, we will work in 
the coulomb gauge and concentrate only on physical transverse gauge field, At{x^ rf) ||Tl],|l3 



Then, the Heisenberg field equation for ^^(x, 77) can be read off from the quadratic part of 
the Lagrangian in the form 



di]' 



■Arix, 7]) - V'^At{x, 7]) + 4c ^(r^) V x At{x, t]) = 0. 



(3.8) 



It is more convenient to decompose the field At{x, rf) into the Fourier mode functions Vy^^{ri) 
in terms of circularly polarized states. 



At{x,vi) 



(Pk 



2{2Tifk 
2{2Tifk 



Mk,ri) 



(3.9) 



where b^^ are destruction operators, and e_j_g are circular polarization unit vectors defined 

Lti( 

+ k''Vik{r])-Akc9{7])Vik{v) = 0, 



in Ref. |TT|. Then the mode equations are 
d'VM 



dr]"^ 

d'V2kiv) 
drf 



+ k''V2k{ri)+Akc9{r])V2k{ri) = ^. 



(3.10) 



with the expectation values given by ||TT| 





" k ' 


j k'^dk coth 


2Ti_ 



d 



^ {\Vik{v)\' - \V2k{v)\' 



(3.11) 
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where we have assumed that at initial time rji, the photons are in local equilibrium at the 
initial temperature Tj. Clearly, the photon mode equations ( p.lO| ) are decoupled in terms of 
the circular polarization mode functions. The axion field acts as the time dependent mass 
term thats triggers photon production. The effective mass terms have opposite signs for 
the two polarizations due to the pseudo-scalar nature of the axion-photon coupling. This 
will lead to producing the different polarized photons in different amounts, resulting in a 
polarization asymmetry in photon emission. The expectation value of the number operator 
for the asymptotic photons with momentum k is given by [|lT 



1 

2k 
1 

1 

4P 



Arik, T]) ■ Ari-k, r]) +PAT{k, r]) ■ Ar^-k, r]) 
k 



- 1 



coth 



coth 



2T, 
k 



2T, 



\VMr + k'\VM\' 



\V2kiv)\^ + k^\V2k{v)\' 



1 

2 + 

1 

2 



N+{k,r]) + N,{k,r]), 



(3.12) 



which is the number of photons with momentum k per unit comoving volume. 



IV. NUMERICAL RESULTS 

In this section, we will compute the photon production from the axionic condensate with 
fa = lO^^GeV, which has a zero-temperature mass rriao — 6.2 x lO^^eV. As we will show 
below, the photon production takes place mainly during the radiation-dominated epoch. So, 
we simply assume a radiation-dominated universe. 

At temperature T, the Hubble parameter is 

(4.1, 

\aj 3 nipi 

where g{T) is the number of effective degrees of freedom at temperature T, rripi is the Planck 
scale, and the cosmic scale factor is 

aiv) = -, (4.2) 

where r/i is the time when the axion field starts to oscillate, defined by a temperature Ti 
such that 3H(Ti) = ma(Ti). As such, r/f^ = ma(Ti)/3. For fa = lO^^GeV, from Eqs. 
and Q we find Ti ~ 0.9 GeV and g{Ti) ~ 60. 

Changing the variable rj into a, the equations of motion ( p.Tj) and (|3.10|) become 



+ _ + 9a2^i-^ sin e - ^ {FF) = , (4.3 

d'^Vi. , de 

d^V2^ .2.. . de 



+ eVa^-t- 4^0—^2? = 0, (4.4) 



da^ da 
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where ^ = krji and 




T » Aqcd, 
T « Aqcd. 



(4.5) 



Note that g{T) does not change significantly from Ti to Aqcd- Henceforth, we approximate 
g{T) ^ g{Ti) ~ 60, where Ti = 0.9 GeV. It is worth to point out that the mode equa- 
tions (|4.4| ) have unstable modes via the spinodal instability for sufficiently low-momentum 
modes with ^ < 4:c\d6/da\, where the effective mass becomes negative. 

To solve Eqs. ( [4.3|) and (|4.4| ), we have to specify the initial conditions for the axion and 
photon fields. The amplitude of the axion field is frozen for r] << rji, i.e., 
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- 
da ' 



as a 



ai « 1. 



For the photon mode functions, we propose 



dVi^ dV, 



2? 



2? 



da 



da 



as a 



a,' << 1. 



(4.6) 



(4.7) 



These initial conditions are physically plausible and simple enough for us to investigate a 
quantitative description of the dynamics. To evaluate the (FF) in Eq. ( p. 11 ) and the photon 
number operator in Eq. ( |3.12| ), we approximate the bose enhancement factor by 



coth 



" k ' 




[el 




= coth 


_2T, 


2r_ 



2r 



(4.8) 



where T = rjiTi > rjiTi ~ 10 , and we are interested in ^ < 100. 

In Fig. ^ we plot the temporal evolution of the axion field and its time derivative by 
choosing a, = 0.01 that corresponds to = 100 Ti. Due to the expansion of the universe, the 
field amplitude decreases with time. However, the rate of change of the amplitude increases 
with time. To understand this, we redefine 6 = 6/a in Eq. ( [4.31) . The (FF) term can be 
neglected, being extremely small as shown in Fig. |^, where we have evaluated the last term 
of Eq. (|4.3| ) denoted by S(a). Then, the equation of motion for 6 when ^ << 1 is given by 



d^9 



2 mliT) ~Q 



0. 



(4.9) 



From Eq. ( [4.5|) , the solutions for 6 are Bessel functions. For T » Aqcd, asymptotically 
e oc cos(0.53a5-^)/a2-35. For T « Aqcd, oc cos(3917 a2)/a0-^ The latter shows that the 
axions behave like non-relativistic matter. 

For the photon production, we calculate the spectral photon number density N{^,a), 
equal to (Nfc(?7)) in Eq. (|3.12| ). It is convenient to define a ratio. 



n(^,a) 



N{^,a)-N{^,a, 



(4.10) 



which is the excess photons above the thermal background. Three snapshots of the ratio 
at a = 0.75, 1.5, and 3 are shown in Fig. 1^. It is interesting to see that the production 
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duration of each non-zero mode is short, and higher- momentum modes are produced at 
later times. This in fact demonstrates a brief exponential growth of the unstable mode due 
to the parametric resonance instability, which is shut-off when the unstable mode has been 
red-shifted out of the unstable band. As a consequence, the photon production is limited 
with an excess photon ratio typically at a level of 10~^. As we have mentioned above, 
the spinodal instability happens for low-momentum photon modes. We demonstrate this 
numerically in Fig. ^ where we have chosen ^ < 0.004 that lie within the spinoidal region 
where ^ < AcldO / da\. The production ratio is also at a level of 10~^, but it is apparent that 
these low-momentum modes are produced during the first oscillating cycle of the axion field. 

The ratio in Fig. |^ can be actually estimated as follows. First, let us find out the 
approximate form for dO/da from Eq. ( |4.9|) and Fig. ||, which is given by 

d6 

— ~ 3.6 a^-^^ cos(0.53 a^'^) for « a < 8.4, (4.11) 
da 

where a = 8.4 is the scale factor when T = Aqcd- At instant a, dO/da is oscillating with an 
effective frequency uj = 0.53 a^'^. Inserting the approximate form ( |4.11|) with uj treated as a 



constant into the mode equation of Vig in Eq. (|4.4|) , and changing variable to 2; = ci;a/2, we 
have 

^ + - 58.1 c-^a'-'' cos{2z)V^^ = 0. (4.12) 



This is the standard Mathieu equations [19|. The widest and most important instability 



is the first parametric resonance that occurs at ^ = uj/2 with a narrow bandwidth 6 ~ 
lA.b ca}'^^ /uj. But actually uj is changing with time. As such, the unstable mode will grow 
exponentially only during a brief period roughly given by 

u5 uj'^5 

^^^A^^^- (4-13) 

Az Aa 

Consequently, this instability leads to the growth of the occupation numbers of the created 
photons by a growth factor, 

g2MA. _ ^ ^ 2/iA;z, (4.14) 

where the growth index fi ~ 6/2. Hence the mode number density is increased by an amount 
approximately given by 

2fiAz ~ 42.6 c^fl-^ ~ 4.6 x lO'^a'^, (4.15) 

where a is the scale factor when the mode Vi^ enters into the resonance band. This estimation 
is of the same order of magnitude as found in the numerical results. Interestingly, the 
dependence of the mode production can also be seen in Fig. |^. Similar estimation can also 
be done for the mode function 

The polarization asymmetry in the produced photons is defined as 

(4.16) 
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We have input the mode solutions for Vi^ and to Eq. ( p.l2|) to calculate S(^,a) ( [4. 161) . 
A plot of the asymmetry versus the momentum at a = 3 is shown in Fig. ^. The numerical 
result shows that the asymmetry is fluctuating about zero as the photon momentum varies. 
The fluctuating amplitude is about 10~^ of the thermal background. Although the asym- 
metry averaged out over a wide range of momenta is nearly zero, at certain momenta the 
produced photons are about 10% circularly polarized. However, subsequently this polariza- 
tion asymmetry will be damped out by photon-electron scatterings in the plasma. 



V. CONCLUSIONS 



First we basically confirm that the photon production via parametric resonance is in- 
validated by the expansion of the universe. Besides, we find that in addition to parametric 
resonance, for long-wavelength photon modes, a new dissipative channel via spinodal insta- 
bility is open. This open channel results in the long-wavelength fluctuations of the photon 
modes. Again, this production is suppressed in the expanding universe. We also observe 
the polarization asymmetry in the produced circularly polarized photons as a result of the 
pseudo-scalar nature of the coupling. However, it will be damped out effectively by the 
plasma. But it is very interesting to see whether it is possible to generate a circular polar- 
ization asymmetry in the production of photons from certain pseudo-scalar fields such that 
it may leave an imprint on the polarization of the cosmic microwave background. As to 
the plasma damping on photon production, we have pointed out the problem in the naive 
approximation that simply introduces the electron plasma frequency to the photon modes. 
We have thus proposed a dynamical and non-equilibrium treatment which should be a bet- 
ter approach to consider the plasma effects. The actual calculations are rather difficult, but 
they certainly deserve further studies. 
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FIGURE CAPTIONS 



Fig. |T[ Time evolution of the axion mean field 9{a) and its time derivative d6{a)/da, where 
a is the cosmic scale factor. 

Fig. 1^. Time evolution of (FF), plotted with the quantity S(a) given by the last term of 
Eq. (13). 

Fig. Ij. Three snapshots of the spectral number density ratio n(^, a), defined in Eq. ( [4.10| ), of 
the photons produced via parametric amplification at a = 0.75, 1.5, and 3. The dimensionless 
quantity ^ = krji, where k is the photon momentum and rji is the conformal time when the 
axion field starts to oscillate. 

Fig. ^. As in Fig. ^, but for low-momentum photons produced via spinodal instability. 

Fig. ^. Circular polarization asymmetry S(^,a), defined in Eq. (|4.16| ), of the photons at 
a = 3 for < ^ < 40. 
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